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ABSTRACT 

We have used the two spectrometers on the Infrared Space Observatory to observe the 
HH 7-11 flow, its red-shifted counterpart, and the candidate exciting source SVS 13, in the 
star formation region NGC 1333. We detect atomic ([Oi]63/xm, [Ol]145/mi, [Siii]34.8/xm, 
[Cn]158/xm) and molecular (H 2 , CO, H 2 0) lines at various positions along the bipolar 
flow. 

Most of the observed lines can be explained in terms of shock-excited emission. In 
particular, our analysis shows that dissociative (J-type) and non-dissociative (C-type) 
shocks are simultaneously present everywhere along both lobes of the flow. We confirm 
the low-excitation nature of the Herbig-Haro nebulosities, with shock velocities v s < 40 — 
50 km s _1 . Toward both lobes of the outflow we find pre-shock densities of no ~ 10 4 cm -3 
for both the J and C components, implying 5o~100uG for Bq(x. n ,' 5 . In the central region 
of the flow, close to the exciting source, the pre-shock density deduced for the C-shock 
component is no ~ 10 5 cm -3 , suggesting a magnetic field ~3 times stronger. We propose 
that the deficiency of gas-phase water in the post C-shock regions is due to freezing onto 
warm grains processed through the J-shock front and traveling along the magnetic field 
lines. The total observed cooling from the dissociative shock components is consistent 
with the power lost by a slow molecular outflow accelerated by a fast neutral Hi wind. 

Finally, the skin of the cloud seen in projection toward the flow appears to be weakly 
photo-ionised by BD +30° 549, the dominant illuminating source of the NGC 1333 reflec- 
tion nebula. 

Subject headings: Stars: formation - (ISM:) Herbig-Haro objects - ISM: individual ob- 
jects: HH7-11 - ISM: molecules - Infrared: ISM: lines and bands 



°Based on observations with ISO, an ESA project with instruments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the United Kingdom) with the participation of ISAS and NASA. 
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1. INTRODUCTION 

Atomic and molecular outflows trace the mass 
loss from protostellar objects, which is a funda- 
mental characteristic of the formation and evolu- 
tion of low mass stars. These outflows are often 
traced by the optical Herbig-Haro (HH) objects, 
shock-excited nebulosities which mark the inter- 
face between outflowing and circumstellar mate- 
rial. One of these systems, which since its discov- 
ery (Herbig 1974 ; Strom et al. 1974) has been 
subjected to a detailed mult i- wavelength analy- 
sis, is HH7-11. The system is relatively bright 
and lies in the very active star forming NGC 1333 
region (Aspin et al. 1994; Bally et al. 1996). 
A distance to the outflow of 350pc (Herbig & 
Jones 1983) is widely adopted in the literature, 
although Cernis (1993) proposes 200pc. Early op- 
tical spectroscopic studies show that the HH7-11 
outflow has a complex velocity field and low exci- 
tation (Solf & B6hm 1987; B6hm & Solf 1990), 
a conclusion further supported by near infrared 
studies (Hartigan et al. 1989; Carr 1993). Early 
on, a well defined CO bipolar outflow was de- 
tected associated with this system (Snell & Ed- 
wards 1981; Bachiller & Cernicharo 1990), which 
also is observed in some other molecules like 
HCO+ and H 2 (Mehringer 1996; Cernicharo et 
al. 1996). 

The HH7-11 outflow has an unusual morphol- 
ogy; the blue outflow lobe is made up of an 
arc-shaped chain of knots, while the red one is 
invisible at optical wavelengths. The red lobe 
is, however, detected in the (l-O)S(l) H 2 line 
at 2.12/zm having a very ragged appearance (see 
Fig 1). Both lobes have a total extension of ~ 
2'. It has been thought that the driving source of 
the outflow is the infrared star SVS 13 (Strom et 
al. 1976), a conclusion partially supported by the 
proper motion measurements of the knots (Herbig 
& Jonesl983) and the source observed outbursts 
(Goodrich 1986; Eisloffel et al. 1991; Liseau et 
al. 1992). The source has a luminosity of ~ 85 Lq 
(Molinari et al. 1993). Recent high angular res- 
olution (~ 0."3) VLA continuum observations at 
3.6 cm suggest another nearby embedded source 
(VLA 3) as a likely candidate, based on its bet- 
ter alignment with the HH string. Interferomet- 
ric observations (Bachiller et al. 1998) at 1.3 and 
3.5 mm with better than 0."2 resolution, however, 
have not confirmed this. The interferometric ob- 
servations have led to the discovery of a second jet 



emanating from a more deeply embedded source 
14."5 away from SVS 13, named SVS 13B (Gross- 
man et al. 1987). 

One of the reasons why the HH7-11 flow 
has been so intensively studied (see e.g. 
Reipurth 1994), is that it was the first system 
showing clear signatures of a high velocity out- 
flow in both neutral and molecular gas tracers 
(Lizano et al. 1988; Rodriguez et al. 1990; Giova- 
nardi et al. 1992). This was a major step forward 
in the interpretation of molecular outflows as be- 
ing driven by faster but more tenuous (than the 
outflow entrained gas) atomic stellar winds (Mas- 
son & Chernin 1993; Raga et al. 1993). This has 
led to a more careful analysis of the energetics and 
shock conditions associated with the ionic/atomic 
and molecular gas outflows (Raga 1991). HH 7-11 
is also one of the few examples where it is possible 
to disentangle the contributions of shock excited 
and fluorescent emission from its near infrared H 2 
spectra (Gredel 1996; Fernandes & Brand 1995; 
Everett 1997). 

In the present study, we take advantage of the 
capabilities of the Infrared Satellite Observatory 
(ISO, Kessler et al. 1996) spectrometers to study 
the mid- and far infrared emission line spectra 
from the HH7-11 red and blue lobes and around 
the driving source SVS 13. The observations are 
described in Sect. 2; the results are presented and 
discussed in Sect. 3 and following, and the main 
conclusions are summarised in Sect. 9. 

2. OBSERVATIONS AND DATA REDUCTION 

We used the two spectrometers on the ISO 
satellite to observe several locations along the 
HH7-11 flow, its optically invisible counterflow 
and the candidate exciting source SVS 13. The 
Long Wavelength Spectrometer (LWS, Clegg et 
al. 1996) was used in its LWS01 grating mode to 
acquire full low resolution (R ~ 200) 43-197/im 
scans with data collected every 1/4 of a resolution 
element (equivalent to ~0.07/im for A < 90/im, 
and to ~0.15/zm for A > 90/im); a total of 19 
scans were collected, corresponding to 38s inte- 
gration time per spectral element. The LWS was 
also used in LWS04 Fabry-Perot (FP) mode to 
collect high (R ~ 8000) resolution scans of the 
[Ol]63/im line; 52 scans sampled at 1/4 of the 
resolution element (i.e., ~ 0.0017 /zm, or 8 km 
s _1 ) were collected, equivalent to an integration 
time of ~ 100s per spectral element. The Short 
Wavelength Spectrometer (SWS, de Graauw et 
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3 26 m 05 s 3 26 m 00 s 3 25 m 55 s 
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Fig. 1. — ISO fields of view in the various observed positions, superimposed onto a continuum- 
subtracted H 2 (l-O)S(l) line map. The dashed circles and solid rectangular apertures represent the 
LWS and SWS fields of view centered on the various observed positions. SVS 13 is at the center of the 
middle LWS (dashed circle) field of view. 



TABLE 1 



Observations 



Object 




a(1950) 


5(1950) 


AOT 


Rev 


Obs. Time 


SVS 13 




03 25 58.1 


+31 05 44.1 


LWS01 


654 


4053 










SWS02 


814 


1730 


HH7 




03 26 02.7 


+31 05 10.2 


LWS01 


654 


3165 










LWS04 


831 


2428 










SWS02 


652 


1464 


HH10 




03 25 59.8 


+31 05 28.8 


SWS02 


847 


1730 


Red Lobe 


1 (RL 1) 


03 25 54.3 


+31 06 34.1 


LWS01 a 


654 


3165 








LWS04 a 


831 


2428 










SWS02 


652 


1464 


Red Lobe 


2 (RL 2) 


03 25 55.4 


+31 06 42.1 


SWS02 


847 


1730 



The field of view also includes the position "Red Lobe 2 
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al. 1996) was used in its SWS02 grating mode 
to observe line scans at medium resolution (R ~ 
2000) for selected wavelength regions covering the 
[Sill]34.8/im and [Nell]12.8/zm fine structure lines, 
and the pure rotational transitions of molecular 
hydrogen from (O-O)S(l) to (0-0) S (7). All rele- 
vant information is summarised in Table 1, in- 
cluding the Astronomical Observation Templates 
used, the revolution number and the total observ- 
ing time of the observations. 

LWS data processed through Off-Line Process- 
ing (OLP), version 7, have been reduced using 
the LWS Interactive Analysis 1 (LIA) Version 7.2. 
The dark current and gain for each detector were 
re-estimated, and the data were recalibrated in 
wavelength, bandpass and flux. The absolute flux 
calibration quoted for LWS in grating mode is 10- 
15% and it is valid for point-like sources since the 
primary calibrator, Uranus, is point-like to the 
LWS beam; however, our sources are not rigor- 
ously point-like and we adopt a more conserva- 
tive number of 20%. Additional processing for 
the FP data (with the LIA routine FP_PROC) 
included a gain correction to compensate for the 
incorrect positioning of the grating during the FP 
observations. The integrated [Ol]63/im line fluxes 
obtained with the FP are about a factor of 2 lower 
than the values measured using the grating. This 
discrepancy is larger than the 30% figure gener- 
ally quoted for the absolute flux calibration of the 
FP (Swinyard et al 1998). We cannot offer any 
explanation for this difference, and in this paper 
we will not use the line fluxes measured with the 
FP. The accuracy of the FP wavelength calibra- 
tion is believed to be better than 1/2 resolution 
element, or ~ 15 km s -1 . 

SWS data were processed using OSIA, the SWS 
Interactive Analysis 2 . Dark currents and photo- 
metric checks were revised; in many cases the for- 
mer were corrupted on a few detectors and were 
re-estimated and subtracted. The March 1998 
bandpass calibration files have been used to pro- 
duce the final spectra. The absolute flux calibra- 
tion for SWS data should be accurate to within 
20%. 

The final steps of data analysis were done us- 
ing the ISO Spectral Analysis Package 3 (ISAP) 
Version 1.5 for both LWS and SWS. Grating 



scans (LWS) and detectors spectra (SWS) were 
averaged using a median clipping algorithm opti- 
mised to flag and discard outliers mainly due to 
transients; line fluxes were estimated by means 
of gaussian fitting (multiple gaussians in case of 
blended lines). The LWS observations toward 
HH 7 and the red lobe (RL 1+2) were heavily 
fringed due to the vicinity of the relatively strong 
continuum source SVS 13; standard techniques 
available under ISAP were used to remove these 
instrumental effects. 

The locations of the observed positions and in- 
strument apertures are shown in Fig. 1 super- 
imposed on an H 2 (1-0)S(1) 2.12/zm continuum- 
subtracted map obtained with the near infrared 
camera at the 60" Mt. Palo mar telescope (Mur- 
phy et al. 1995). The apertures of the two ISO 
spectrometers are quite different. The SWS focal 
plane aperture is 14"x20" for all of the detected 
lines except for [Sin]34.8/_m (20"x33") so that, 
at each pointing, the contamination from nearby 
knots should be negligible. The LWS aperture is 
rather large (0 ~80", Swinyard et al. 1998); the 
pointing on SVS 13 also includes HH11 and 10 
and two sources VLA 3 and SVS 13B (see Sect. 1), 
while that centered on HH 7 also includes HH8, 9 
and 10. 

3. RESULTS 

Table 2 presents the fluxes for the lines detected 
toward each location. One sigma upper limits 
are given only for lines observed with dedicated 
SWS02 "line scans" AOTs but not detected (see 
Tab. 1); a horizontal dash means that no obser- 
vation is available for that particular line. 

Detected lines are also plotted in Figs. 2a, b, 
c and d; most of them are detected everywhere. 
Exceptions are [Sill]34.8/im, which is not detected 
toward RL2, ortho-B. 2 which is not detected to- 
ward RL1+2, and [Nell]12.8;um which is not de- 
tected anywhere. All lines, except [Sill]34.8/im 
and the H 2 lines, are stronger on the SVS 13 po- 
sition. 

The FP spectra of the [Ol]63/im line toward 
HH 7 and RL 1+2 are presented in Fig. 3. The 
lines are resolved (Avfp ~30 km s~ x ) with de- 
convolved FWHM of ~50 km s _1 and ~25 km s _1 
for the two positions; these widths are consistent 
with the velocity field traced by the [S II]6717A 

■"■LIA is available at http://www.ipac.caltech.edu/iso/lws/lia/lia.html 

2 OSIA is available at http://www.mpe.mpg.de/www_ir/ISO/observer/osia/osia.html 

3 ISAP is available at http://www.ipac.caltech.edu/iso/isap/isap.html 
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TABLE 2 



Observed Line Fluxes a 



Line 


SVS 13 


HH 7 


HH 10 


Red Lobe 1 


Red Lobe 2 








nes 






(0-0)S(l) d 


4.0(0.5) 


2.6(0.5) 


5.5(0.5) 


4.5(0.7) 


4.2(0.5) 


(0-0)S(2) c 


7.9(2.6) 


13.3(1.9) 


6.9(1.2) 


7.6(2.1) 


8.3(2.5) 


(0-0)S(3) c 


9.3(0.8) 


10.2(0.7) 


12.5(1.1) 


10.6(0.9) 


7.8(0.8) 


(0-0)S(4) c 


4.2(0.6) 


10.9(1.0) 


5.3(0.6) 


4.2(0.8) 


3.1(0.7) 


(0-0)S(5) c 


9.9(2.0) 


15.5(1.9) 


9.9(1.5) 


8.7(2.0) 


11.9(2.4) 


(0-0)S(6) c 


<2 




<3 




<3 


(0-0)S(7) c 


<4 




<5 




<5 


[Nell]12.8/zm d 


<0.7 


<0.9 


<0.7 


<0.9 


<0.8 


[Siii]34.8/xms 


7.4(2.3) 


3.6(1.1) 


6.8(2.2) 


5.5(1.0) 


<3 



LWS Lines 

CO 20-19* 5.4(1. 7) b = 
CO 19-18 f 13.7(2.0) b 
CO 18-17 f 13.0(4.6) b 

CO 17-16 f 10.1(1.7) b 8.4(3.0) c - 6.6(1.2) 

CO 16-15 f 10.1(1.7) b 5.7(1.5)° - 3.9(0.6) 

CO 15-14 f 13.2(6.3) b 8.5(2.8) c 

CO 14-13 f 10.1(2.7) b 6.6(2. l) c 

o-H 2 3(0,3)-2(l,2) f 12.2(6.3) b 6.5(2.8) c 

o-H 2 2(l,2)-l(0,l) f 6.8(3.8) b 7.9(1.5) c 

[Ol]63/xm f 180(2) b 131(3) c - 109(3) 

[Ol]145/zm f 8.0(4.0) b 4.6(1.0) c 7.9(1.3) 

[Cll]158^m f 21.3(1.4) b 12.5(1. l) c - 15.2(0.7) 

a Units of 10~ 20 W cm _2 and la uncertainties in parenthesis; upper limits are at the la level. 

A horizontal dash means that no observation is available for that particular line. 
h Includes HH 10 and HH11. 

c Includes HH 8 and; HH 9 and HH 10 are at the edge of the LWS beam. 
d Focal plane aperture is 14" x 27". 
e Focal plane aperture is 14" x 20". 
/ Focal plane aperture radius is 40". 
9 Focal plane aperture is 20" x 33". 



HH 7A - [01] S5V15 - [Ol] RL1+2 - [01] 




^ 61 62 63 64 65 61 62 63 64 65 61 62 63 64 65 

Wavelength (fJ.rn) 

Fig. 2. — a. [Ol]63/im lines detected in the LWS full grating scans on the three observed positions 
(dashed circles in Fig. 1). 
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HH 7A 




E 4 - 
E 3 - 






150 160 

Wavelength (/im) 



Fig. 2. — b. Lines detected between 130/im and 190/im in the LWS full grating scans on the three 
observed positions (dashed circles in Fig. 1). 



HH7A - S(1) 



16.98 17.03 17.1 

HH10 - S(1) 



16.98 17.03 17.08 
SSV13 - S(1) 



16.98 17.03 17.08 
RL1 - S(1) 



16.98 17.03 17.08 
RL2 - S(1) 



16.98 17.03 17.08 



HH7A - S(2) 



HH7A - S(3) 





12.23 12.28 12.33 9.63 9.66 9.69 

HH10 - S(2) HH10 - S(3) 



^^^^^^ 




12.73 12.28 17.33 9.5.5 9.56 9.69 

SSV13 - S(2) SSV13 - S(3) 




12.23 12.28 12.33 

RL1 - S(2) 



9.63 9.66 9.69 
RL1 - S(3) 




12.23 12.28 12.33 

RL2 - S(2) 



9.63 9.66 9.69 
RL2 - S(3) 



17.73 17.78 17.33 9.63 9.66 9.69 

Wavelength (/im) 



HH7A - S(4) 



8.00 8.03 8.06 

HH10 - S(4) 



8.00 8.03 8.06 
SSV13 - S(4) 



8.00 8.03 8.06 

RL1 - S(4) 



8.00 8.03 8.06 

RL2 - S(4) 



^^^^^^^ 



8.00 8.03 8.06 



HH7A - S(5) 




6.90 6.97 6.94 
HH10 - S(5) 




6.90 6.92 6.94 
SSV13 - S(5) 




6.90 6.92 6.94 
RL1 - S(5) 



6.90 6.97 6.94 

RL2 - S(5) 




6.90 6.92 6.94 



Fig. 2. — c. H 2 pure rotational lines detected in the SWS line grating scans on the five observed 
positions (full rectangles in Fig. 1). 
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HH 7A - [Sill] HH 10 - [Sill] SSV1 3 - [Sill] RL1 - [Sill] 




34.70 34.80 34.90 34.70 34.80 34.90 34.70 34.80 34.90 34.70 34.80 34.90 

Wavelength (fj.m) 



Fig. 2. — d. [Siii]34.8//m lines detected in the SWS line grating scans on four out of five observed 
positions (full rectangles in Fig. 1). 



1x10 
8X1CT 




-2xio 1 7 r . ' .. . . i . . . . i . . . . i . . ' . ; . 

-100 -50 50 100 

Velocity (km s~ ) 



Fig. 3.— [Ol]63/xm line detected with the LWS Fabry-Perot toward HH7 (full line) and RL1+RL2 
(dashed line); the dotted vertical line indicates the cloud systemic velocity, 8.3 km s -1 . 



line (Stapelfeldt 1991) toward HH7, and suggest 
that [Ol]63/im originates from the flow material. 

4. THE GAS PHYSICAL PARAMETERS 
4.1. H 2 

The observed H 2 lines arise from quadrupole 
rotational transitions of the ground vibrational 
state so that they are likely to be optically thin. 
The radiative lifetimes range from 0.54 yrs for 
level v—0, J=7 where line (0-0) S (5) originates, to 
~ 1000 yrs for level v—0, J=2 which is the up- 
per state for (0-0) S(0) transition; thus these lev- 
els are collisionally populated, and a simple LTE 
analysis is adequate to interpret the data (e.g., 
Gredel 1994). Einstein coefficients and wavenum- 
bers were taken from Black & Dalgarno (1976) 
and Dabrowski (1984); ortho/para=3 has been 
assumed. We dereddened the H 2 line fluxes us- 
ing the visual extinction given by Gredel (1996) 
for the individual HH knots, and the Rieke & 
Lebofsky (1985) extinction curve. Dereddened 
line fluxes were used to produce the Boltzmann 
plots shown in Fig. 4; a single temperature fit is 

4 We believe that Table 3 of Gredel (1996) incorrectly r< 
note c of that Table, the total column densities should be a 
(R. Gredel, priv. comm.) the IRSPEC slit was not aligned 



in good agreement with the observed fluxes at all 
positions. The solid angles of line emitting re- 
gions are arbitrarily set to the equivalent SWS 
focal plane aperture, i.e. 14" x 20" or 6.6 xlO -9 
sr (valid for all detected H 2 lines but S(l), where 
the solid angle is ~30% higher), assuming a beam 
filling factor of 1. 

The H 2 temperatures and the SWS beam av- 
eraged column densities, together with la uncer- 
tainties are listed in Table 3. The H 2 column den- 
sities should be considered as lower limits because 
the beam filling factor can be less than 1. The H 2 
temperatures vary by less than 15% along the flow 
and appear to trace an H 2 component (which we 
call "warm") which is colder than the one iden- 
tified by Gredel (1996) via higher excitation H 2 
ro-vibrational lines (2100 K< T <2750 K, which 
we will call "hot"). The column densities for the 
warm H 2 (see Table 3, col. 3) are on average two 
orders of magnitude higher than those of the hot 
H 2 (Gredel 1996) and the difference cannot be ex- 
plained by the smaller (a factor 2) emitting areas 
adopted by Gredel 4 . This suggests that the two 

ports the H2 column densities; based on the explanation in 
actor ten higher than reported in its column 7. Furthermore 
exactly on the peaks of the individual knots, suggesting that 
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1000 2000 3000 4000 5000 6000 7000 

E(J up ) (K) 

Fig. 4. — H 2 excitation diagrams for the five observed positions. In each panel, the abscissa represents 
the energy of the level in K, while Nj/gj is in ordinata; offsets of -10, -5, 5 and 10 dex have been applied 
to data of SVS 13, HH 10, RL 1 and RL 2 while no offset has been applied to the HH 7 data. 



TABLE 3 



H 2 Physical Parameters 



Object 


Th 2 


N(H 2 ) 




(K) 


(10 20 cm" 2 ) 


SVS 13 


560(20) 


0.7(0.1) 


HH7 


670(20) 


0.5(0.1) 


HH10 


540(10) 


1.0(0.1) 


RL 1 


540(20) 


0.9(0.2) 


RL 2 


560(20) 


0.8(0.1) 
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H 2 components may be actually distinct. 
4.2. CO 

CO rotational lines have been detected at all 
three positions observed with LWS, although the 
intensity and number of detected lines are high- 
est toward SVS 13. More CO lines than those 
reported in Table 2 are marginally visible in our 
spectra toward the two flow positions, but it is im- 
possible to reliably determine their flux because of 
the heavily fringed LWS spectra. Contrary to the 
H 2 lines, the CO rotational spectrum arises from 
dipole transitions and the simplifying assump- 
tions of optically thin lines cannot be adopted a 
priori. We analysed the CO lines using the LVG 
model described in Ceccarelli et al. (1998). Un- 
der NLTE conditions the line ratios depend both 
on the gas temperature and density, as well as on 
the CO column density if the lines are optically 
thick. The observed CO lines cannot constrain 
the three parameters simultaneously. The distri- 
bution of CO line fluxes vs J toward SVS 13 is 
essentially flat, and we find a range of physical 
conditions which are consistent with our obser- 
vations. We find, as extreme cases, a 'cold' so- 
lution with T~350 K and n ~ 10 6 cm -3 , and a 
'warm' solution with T~900 K and n ~ 10 5 cm -3 ; 
the model also predicts optically thin CO lines. 
Interestingly, this temperature range is centered 
around the value independently derived from H 2 
lines, which suggests that the CO and H 2 emis- 
sion come from the same region. For T=560 K 
(the H 2 temperature toward SVS 13, see Table 3) 
we obtain n ~ 4 x 10 5 cm" 3 . The CO line ra- 
tios toward HH7 and RL 1+2 imply temperatures 
and densities similar to those around SVS 13, al- 
though the lower absolute line fluxes would sug- 
gest column densities or beam filling factor about 
a factor two lower. Using the SWS beam solid 
angle we derive N(CO)~ 2 x 10 16 cm -2 toward 
the SVS 13 LWS pointing, giving a direct mea- 
surement of the CO abundance in this warm gas 
of [CO]/[H]~ 1.2 x 10" 4 . Based on this number 
it seems that CO accounts for essentially all of 
the gas-phase carbon in the interstellar medium 
(Cardelli et al. 1996). Recently Lefloch et al. 
(1998) found evidence of CO depletion toward the 
core of SVS 13, likely due to condensation onto 
grains in the dense environments close to SVS 13. 
We find no evidence of such a depletion in the 
high-J CO line emitting region, suggesting that 

quoted column densities should more conservatively be considered as lower limits. 



all of the CO locked into the grain mantles has 
been returned to the gas-phase. 

4.3. H 2 

As far as H 2 is concerned, since only the 
3(0,3)-2(l,2) 174.6/im and 2(1,2)-1(0,1) 179.5/zm 
lines of ortho-R 2 have been detected in SVS 13 
and HH 7, it is clear that water plays only a sec- 
ondary role in gas cooling (see Table 4) compared 
to some other cases of outflow exciting sources 
such as IC 1396N (Molinari et al. 1998) or L1448- 
mm (Nisini et al. 1999). This result goes against 
earlier suggestions based on previous A=1.67 mm 
water line observations (Cernicharo et al. 1996), 
that the emission originates in dense (10 6 cm" 3 ) 
shocked material with a water abundance com- 
parable to that of CO. It should be noted, how- 
ever, that Cernicharo et al. assumed a gas tem- 
perature of 50 K, while our observations of H 2 
and CO lines clearly indicate that the tempera- 
ture of the molecular material is a factor of 10 
higher. The maps presented by Cernicharo et al. 
show that water emission is concentrated within 
~ 10" region centered on HH 11. Adopting T = 
560 K and n = 4xl0 5 cm -3 , compatible with the 
conditions for H 2 and CO, our model fit (Cecca- 
relli et al. 1998) to ISO and millimeter water lines 
predicts optically thick lines and a water column 
density N(H 2 0)~ 10 15 cm" 2 , about a factor 20 
lower than that estimated by Cernicharo et al. 
An abundance [H 2 0]/[H 2 ] <6x 10" 6 is also de- 
rived assuming that the lines from two molecules 
come from the same gas component; we will prove 
this assumption in Sect. 6. 

5. THE PHOTO-DISSOCIATION REGION 

In spite of the shock-excited nature of the HH 
objects, it is likely that a non-negligible con- 
tribution to the observed line emission actually 
come from an extended Photo-Dissociation Re- 
gion (PDR, Tielens & Hollenbach 1985) compo- 
nent associated with the NGC 1333 cloud. This 
is suggested by low dispersion ISO-LWS observa- 
tions (Caux et al. 2000) at various positions in 
the NGC 1333 cloud, which show a widespread 
[Ol]63;um and [Cn]158;um emission; this extended 
component seems to account for ~20% of the 
[Ol]63/im, and all of the [Cn]158;um emission 
we see from our LWS pointings. Our FP data 
(Sect. 3) confirm that only a small fraction of the 
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be as high as 10 5 K, resulting in complete molec- 
ular dissociation. In C-shocks, the temperature 
rarely exceeds a few thousand degrees and molec- 
ular material can survive. These two very differ- 
ent physical scenarios produce distinctive signa- 
tures in terms of cooling ratios between different 
species, and of line ratios within the same species. 
We will show that our data of the HH 7-11 region 
depict a complex situation where the two types 
of shock coexist; to help the discussion below, we 
report in Table 4 the total cooling rates in the 
various species as, when applicable, derived from 
the models used to estimate their physical param- 
eters. 

6.1. J-shocks 

First of all we note that the non-detection of the 
[Nell]12.8/iin line down to a level of 10~ 20 W cm~ 2 
confirms the low-excitation nature of the HH 7-11 
chain and suggests (Hollenbach & McKee 1989, 
hereafter HM89) a shock velocity v s < 40 — 50 
km s -1 , depending of the pre-shock density, in 
excellent agreement with our high-resolution FP 
[Ol]63/im observations (Sect. 3) and with esti- 
mates from optical spectroscopy (Bohm, Brugel 
& Olmsted 1983; Solf & B6hm 1987). This up- 
per limit on the shock velocity excludes strong 
shocks, yet the detection of [Sill]34.8/im requires 
the presence of a dissociative shock component, 
since negligible ionization is expected from a C- 
shock (HM89). Such a component would also ex- 
plain the observed [Ol]63;um cooling, since the 
latter is expected to be the main coolant in J- 
shocks 7 . In order to compare the [Ol]63/im with 
the [Siii]34.8//m line for HH 7 and SVS 13, we need 
to determine the fraction of [Ol]63;um emission 
due to other HH objects falling within the LWS 
beam. 

The LWS beam centered on HH 7 also contains 
HH8, 9 and 10; since the excitation conditions 
for the different HH objects along the flow do 
not show dramatic variations (Hartigan, Curiel 
& Raymond 1989), we choose to use the de- 
tected [Sill]34.8//m lines toward HH 7 and HH10 

5 Thc "PDR Toolbox" is available at http://dustem.astro.umd.edu and contains downloadable FIR lines diagnostic 
information about PDRs. The tool has been created by L. Mundy, M. Wolfire, S. Lord and M. Pound, and it is based on 
the new PDR models of Kaufman et al. (1999). 

6 The ISO archive is available at http://pma.iso.vilspa.esa.es 

7 Copious [Oi]63/im can also be produced in non-dissociative shocks. In this case however, the presence of H2 would 
allow the incorporation of O into water, via the chain of endothermic reactions O + H2 — > OH + H and OH + H2 — > H2O 
+ H. This chain has an activation energy of T^ 220 K, lower than the H 2 temperature estimated from the rotational lines 
(between 540 and 670 K, see Table 3); we conclude that the observed [Oi]63^m emission cannot originate in C-shocks. 



[Ol]63yum line may come from a quiescent PDR 
component. 

Using the PDR Toolboafi, the observed 
[Cii]158/xm emission requires a relatively faint 
FUV irradiation level of G ^$ 10 in units of lo- 
cal Galactic FUV flux (Habing 1968). Although 
SVS 13 might be able to provide the required 
FUV field, the widespread [Cll]158//m emission 
seen by Caux et al. (2000) clearly suggest an 
external irradiation source. A natural candidate 
is BD +30° 549, the B6 star responsible for the il- 
lumination of the NGC 1333 reflection nebulosity 
(Harvey et al. 1984); located ~ 0.8 pc N-NE of the 
HH7-11 area, it can certainly provide the needed 
G < 10 FUV field. In this regime, the PDR 
surface temperature does not exceed ~ 100 K 
(Kaufman et al. 1999, Timmermann et al. 1996, 
Kemper et al. 1999, Liseau et al. 1999) which 
excludes a PDR origin for the [Sill]34.8/im line 
(Hollenbach, Takahashi & Tielens 1991), the H2 
and the CO lines. Indeed, the temperature of 
the molecular material is at least 5 times higher, 
and the molecular emission does not appear to be 
extended as one would expect for a PDR origin; 
ISOCAM-CVF near-IR imaging spectroscopy of 
the HH7-11 region extracted from the public ISO 
data archive 6 shows that the emission from the 
same H 2 rotational lines observed with SWS is 
concentrated along the flow and peaks in corre- 
spondence of visible HH objects (Noriega-Crespo 
et al. 2000), obviously favouring a shock origin 
for these lines. 

6. SHOCKS ALONG THE FLOW 

It is known that the nature of the shock exci- 
tation is dramatically influenced by the presence 
of a magnetic field component perpendicular to 
the shock velocity which prepares the up-stream 
medium and smoothes out the effect of the front 
passage. The differences in the physical condi- 
tions of shocked e; clS clXG such that two distinct 
classes of shocks, C(ontinuous) and J(ump) have 
been idealised (Draine 1980). In a J-shock the 
temperature reached by the shocked material de- 
pends on the square of the shock velocity and can 
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TABLE 4 
Cooling rates' 



Species 


HH7 


HH10 


SVS13 


RL 1 RL 2 


[OlJ63/im 
CO 


4.8 
2.0 




6.6 
3.0 


4.0 
1.4 


H 2 


0.7 




0.7 




H 2 

[Sill]34.8um 


2.3 


1.9 


1.6 


1.7 1.6 


0.13 


0.25 


0.27 


0.20 



Units of 1CT 2 L, 



as weights to estimate the [Ol]63um cooling in- 
trinsic to HH 7. We conservatively assign to HH 8 
and HH9 the same [Sill]34.8/xm flux measured to- 
ward HH 7, and we also consider that the weight- 
ing of the LWS beam profile decreases the con- 
tribution of HH8-9 and HH 10 to the [Ol]63um 
measured on the HH 7 pointing by 10% and 50%, 
respectively. Similar estimates can be done for 
the SVS 13 pointing, where again HH 10 con- 
tributes at a 50% level. The HHll's contribu- 
tion does not need to be disentangled since it 
also contributes to the [Siii]34.8//m line of SVS 13. 
No correction is required for the RL 1+2 point- 
ing. Taking the above into account, we obtain 
[Ol]63um/[Sill]34.8um line ratios between 15 and 
20, suggesting a pre-shock density no ~ 10 4 cm -3 
(HM89). The models can also reproduce the ab- 
solute fluxes, provided that the emission solid an- 
gles are ~ 5" in diameter with filling factor of 1; 
more intense Ol and Sill lines on the central po- 
sition are likely due to slightly larger solid angles 
and/or filling factors. 

As concerns molecular cooling the main con- 
tribution to H 2 emission in J-shocks is predicted 
(HM89) to come from material excited by FUV 
or H2-formation pumping. Indeed, Fernandes & 
Brand (1995) propose a fluorescent origin for the 
near-IR H 2 lines toward HH 7. For the (0-0) S lines 
however, the predicted line ratios are not repro- 
duced by the observations. In particular the S(7) 
line, which is predicted to be always brighter than 
the S(5) line irrespectively of pre-shock density 
and shock velocity, is not detected at all in our 
SWS spectra 8 . 



6.2. C-shocks 

When interpreted in terms of non-dissociative 
shock, the pure rotational H 2 lines provides a sen- 
sitive probe for the shock velocity, given the four 
orders of magnitudes of dynamical range spanned 
by their line ratios. This is shown in Fig. 5, which 
presents the observed H 2 line ratios superimposed 
on a grid of C-shock models from Kaufman & 
Neufeld (1996). 

Our observations seem to place quite stringent 
boundaries to the shock velocity, 15 < v s < 20 
km s _1 . The corresponding gas temperature, 
~ 600 K (Kaufman & Neufeld 1996), is in good 
agreement with the values derived using simple 
LTE analysis (Sect. 4), as expected since the 
molecular gas is collisionally heated by the C- 
shock front. Observed absolute fluxes can also 
be reproduced by the model as long as the emis- 
sion solid angle does not exceed few arcsecs in 
diameter. We note that such a low velocity shock 
would also produce very faint S(6) and S(7) lines. 
Once the shock velocity is determined, we can use 
the CO/H 2 cooling ratio to estimate the pre-shock 
density. Fig. 6 presents such a diagnostic dia- 
gram. 

The CO line fluxes in the three LWS point- 
ings should be corrected for contamination by 
the other HH objects. In similar way as for the 
[Ol]63um/[Sill]34.8um ratio, we use the H 2 cool- 
ings in different positions as weights to split the 
individual contributions to the CO cooling. The 
corrected cooling ratios are reported as dashed 
lines in Fig. 6. We see that a shock velocity 
15 < v s < 20 km s _1 would correspond to a pre- 
shock densities of n ~ 10 4 cm" 3 for HH7 and 



8 The above mentioned ISOCAM-CVF observations (at lower spectral resolution than the SWS) actually detected the 
S(7) but with a flux about 4 times lower than the S(5) line flux, compatibly with our SWS upper limits on the S(7) line 
(Noriega-Crespo et al. 2000). 
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Fig. 5. — C-shock diagnostic ratios for the H 2 pure rotational lines. The solid and dashed lines repre- 
sent iso-density and iso-velocity levels respectively. The levels are labeled in cm" 3 (powers of ten) and 
in km s -1 . The crosses represent the H 2 line ratios for the five positions observed with the SWS. 




Fig. 6. — CO/H2 cooling ratio as predicted for C-shock models (Kaufman & Neufeld 1996) as a 
function of pre-shock density and shock velocity. The dashed lines indicate the values observed in the 
different LWS pointings. 
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Fig. 7. — Decimal logarithm of the CO/H2O cooling ratio as predicted for C-shock models (Kaufman 
& Neufeld 1996) as a function of pre-shock density and shock velocity. The dashed lines indicate the 
values observed toward HH7 and SVS 13 (water emission was not detected toward RL 1+2). 



RL 1+2, and n ~ 10 5 cm^ 3 for SVS 13. 

Water is an important diagnostic for C-shock 
models. H 2 cooling is a fast function of the 
shock velocity since free atomic oxygen in the 
post-shock gas is expected to be incorporated 
into water as soon as the temperature rises above 
~ 200 K, roughly corresponding to v s ~ 10 km s _1 
(Draine, Roberge & Dalgarno 1983, Kaufman & 
Neufeld 1996). In our case then, where the tem- 
perature is ~600 K, cooling via H2O rotational 
transitions is expected to be dominant with re- 
spect to that of CO. Instead, the diagnostic dia- 
gram in Fig. 7 shows that the observed CO/H 2 
cooling ratio is consistent with shock velocities 
v s < 12 km s _1 , which cannot justify the observed 
H 2 temperatures. This suggests that most of the 
expected gas phase water is missing. We will pro- 
pose an explanation for this result in Sect. 7.1 

7. THE SHOCKS AND THE HERBIG-HARO 
OBJECTS 

We have shown that HH7-11 is a complex re- 
gion where different line emission mechanisms are 
simultaneously at work. Notwithstanding the 
poor spatial resolution of our data, the variety of 
different spectral signatures detected allow us to 



draw the following physical scenario of the region. 
7.1. HH7 

Starting from HH 7, the C-shock conditions di- 
agnosed by the molecular emission clearly re- 
quire the presence of a magnetic field B +v sc 
(Draine 1980), where v Sc is the shock velocity of 
the C component. Assuming the standard scaling 
law (Draine 1980) for the value of B 

B ~ b n° 5 /j,G (1) 

where b ~ 1 in the interstellar medium, the es- 
timated n ~ 10 4 cm~ 3 pre-shock density for the 
C-shock component on HH7 gives B ~ 100/zG. 
Variations of a factor 3 in each direction are nev- 
ertheless possible, since b can vary between 0.3 
and 3 in molecular clouds (HM89). For similar 
B and no, Draine, Roberge & Dalgarno (1983) 
have shown that the C^J transition occurs at 
v s ~50 km s _1 , similar to the upper limit set on v Sj 
(the J-shock velocity) by our observations. This 
means that a significant Bo component transverse 
to v Sj would smooth the shock front to a C type 
or, equivalently, that our observed J-shock com- 
ponent can only exist as long as v S/ ||B . When 
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the bowshock-like morphology of HH 7 is also con- 
sidered, as revealed by HST NICMOS 2.12um im- 
ages (A. Cotera, priv. comm.), then a very sim- 
ple scenario emerges where HH 7 is immersed in 
a B roughly parallel to the flow axis. At the tip 
of the bow v s ~||B and J-shock conditions are 
present; [Sill]34.8um and most of [Ol]63/xm flux 
arise from this region, and favourable conditions 
also exist to originate the "hot", T~ 2100 K, H 2 
emission detected by Gredel (1996). Along the 
sides of the bow, v s becomes nearly perpendicular 
to B , creating favourable conditions for C-type 
shocks. This picture also provides a plausible ex- 
planation to the problem of the missing water: 
the water is actually produced in the C-shocks, 
but rapidly condenses onto dust grains and dis- 
appears from the gas-phase. This possibility is 
suggested by the recent discovery with ISO-LWS 
(Molinari et al. 1999) of crystalline water ice to- 
ward HH7. The deduced water abundance (in 
solid state form) is comparable to the interstel- 
lar oxygen abundance, which is what the mod- 
els would predict for gas-phase water in C-shocks; 
furthermore, the fact that the ice is in crystalline 
form requires grain temperatures of the order 
of 100 K, only attainable in dissociative shocks 
(HM89; Draine, Roberge & Dalgarno 1983). Af- 
ter being heated, the grains would be efficiently 
transported in the post C-shock regions along the 
B lines, directed parallel to the flow. For the gas 
parameters we derived, the grains are efficiently 
coupled to the magnetic field and are not signifi- 
cantly decelerated by collisions with the neutrals 
(Draine, Roberge & Dalgarno 1983). 

7.2. HH10 

From the viewpoint of the line emission proper- 
ties, HH 10 and HH 7 appear very similar objects. 
Although no LWS data were specifically collected 
toward HH 10, the H 2 and the [Sill]34.8um lines 
still argue in favour of a dual J+C shock nature. 
From the morphological point of view, HH 10 ap- 
pears as an irregular blob in the 2.12/xm image of 
Fig. 1. Higher spatial resolution images in Ha and 
[Sll]6717+3lA from HST (unpublished archival 
data) resolved HH 10 into a double filamentary 
structure whose N-S orientation does not appear 
to be related to the flow direction. However, the 
Ha/[Sll]6717+3lA emission ratio is higher in the 
NW part of this HH object, facing toward SVS 13. 
A higher ratio implies higher excitation (stronger 
shocks) conditions, which are likely to be traced 



by the [Siii]34.8//m line and by the 'hot', T~ 2200 
K, H 2 component (Gredel 1996). If we assume 
that Bo maintains its direction downstream along 
the flow, as suggested for HH7 (see above), then 
we speculate that the morphology of HH 10 could 
correspond to irregularities or corrugations in the 
shocked walls of the cavity which is excavated by 
the flow. Such irregularities might originate from 
instabilities at the flow-cavity interface, as pro- 
posed by Liseau, Sandell & Knee (1988), although 
there are no detailed numerical simulations of the 
process. 

7.3. HH 11 and SVS 13 

The interpretation for HH11 and SVS 13 is 
more complicated because both sources are con- 
tained in the fields of view of the SWS and LWS 
instruments. Pre-shock densities no ~ 10 5 cm" 3 
are here found for the C component; this is a fac- 
tor ten higher than in the other positions of the 
flow, which is not surprising given the close prox- 
imity of the origin of the flow. In these conditions 
the magnitude of the magnetic field (Eq. 1) could 
range from ~ 300uG to ~lmG. Such high B val- 
ues have been claimed by Hartigan, Curiel & Ray- 
mond (1989) to justify the relative faintness of 
HH 11 in 2.12um H 2 images (see also Fig. 1). The 
possibility that FIR line emission toward SVS 13 
arises in a collapsing envelope around is not rele- 
vant here because the predicted [Ol]63um line flux 
(Ceccarelli, Hollenbach & Tielens 1996) is about 
30 times lower than actually observed, while CO 
and H 2 lines are below the detection limit for the 
present observations. Finally, there is the possi- 
bility that a fraction of the line fluxes measured 
with the LWS originates from the recently discov- 
ered embedded outflow source SVS 13B (Bachiller 
et al. 1998). 

7.4. The Red- shifted Lobe 

[Ol]63um, [Ol]145um and [Sill]34.8/xm, to- 
gether with a complement of H 2 and CO lines, 
have been detected toward the receding lobe, and 
define shock conditions which are not dramati- 
cally different from those present on the blue lobe. 
It is well known that no optical emission is de- 
tected toward the red lobe, and also published 
images in the H 2 2.12//m line (Fig. 1, Garden 
et al. 1990, Hodapp & Ladd 1995) clearly show 
fainter emission there. Higher values of dust ex- 
tinction with respect to the blue lobe have been 
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invoked as an explanation for this asymmetry. 
Our observations, which trace similar shock con- 
ditions for the two lobes, tend to support this pos- 
sibility. 

8. THE SHOCKS AND THE MOLECULAR 
OUTFLOW 

If J-shocks between stellar winds and ambi- 
ent material are responsible for the acceleration 
of the molecular outflow, a correlation (Hollen- 
bach 1985) is expected between the outflow mass 
loss rate and the flux of the [Ol]63/xm line, which 
is the dominant coolant in such shocks. The 
predicted mass loss rate, based on the observed 
[Ol]63/im cooling from HH7, is 4.8 x 1CT 6 M 
yr" 1 (see also Cohen et al. 1988, Ceccarelli et 
al. 1997), in good agreement with the mass loss 
rate estimated by Lizano et al. (1988, see also Ro- 
driguez et al. 1990) for the fast Hi wind believed 
to be responsible for the acceleration of the slow 
CO outflow (Snell & Edwards 1981) . This fast 
neutral wind was also confirmed with CO obser- 
vations by Bachiller & Cernicharo (1990). 

Assuming momentum balance at the interface 
between the the wind and the ambient medium, 
Davis & Eisloffel (1995, 1996) derived a simple 
relationship between the mechanical power of the 
wind L w and the power radiated by the shock 

Lrad- 



J rad 



1-^ 

V m . 



(2) 



Since the working surfaces where the winds 
impact the medium are traced by the J-shocks, 
we make the assumption L ra d ~ Loi + Lsm + 
L H2 _ rovib . The total cooling due to the near-IR 
H 2 vibrational lines measured by Gredel (1996) 
along the flow is ~ 5 x 10~ 3 L Q ; we assume a 
slit width of few arcsecs, so we will conservatively 
multiply the observed value by 10 to allow for the 
extension of the HH objects. We estimate the 
mechanical power of the Hi wind according to: 



l M w v 2 w 

2 Tdyn 



(3) 



Using the Hi parameters from Lizano et al. 
(1988), we obtain L w ~ 1.9 L & . Eq. (2) then pro- 
vides v s /v w ~ 0.6, or f S/ ~ 36 km s" 1 for an aver- 
age Hi wind velocity of v w ~60 km s" 1 (Lizano et 
al. 1988), in excellent agreement with our results. 



9. SUMMARY 

The HH7-11 flow, together with its red-shifted 
counterpart and SVS 13 (the candidate excit- 
ing source) have been studied via atomic, ionic 
and molecular spectroscopy. A complex scenario 
emerges, where: 

1. we have detected atomic ([Ol]63/im, 
[Ol]145/im), ionic ([Cll]158/im, [Sill]34.8//m) 
and molecular (H 2 , CO and H 2 0) lines 
along the flow (both lobes) and toward 
SVS 13. 

2. the low-excitation shock nature of the HH 
nebulosities along the flow is confirmed. 
Spectral signatures of C and J shocks are 
ubiquitously found along the HH7-11 flow 
and its red-shifted counterpart. Our esti- 
mates for the shock velocities are v Sj < 40 — 
50 km s" 1 and 15 < v Sc < 20 km s _1 . The 
pre-shock density is ~ 10 4 cm -3 toward the 
blue and the red lobe; for the C component 
only, we find Hq ~ 10 5 cm" 3 at the location 
of SVS 13. 

3. there is indirect evidence for an ordered B 
field oriented parallel to the direction of the 
flow. The magnitude of the magnetic field 
is £? ~ lOOyuG on the lobes, increasing to 
~ 300/iG at the position of the flow origin; 
these figures, however, can vary of a factor 
3 in each direction. 

4. the gas-phase in the post C-shock region 
is deficient in H 2 0. We presented evidence 
that this may be due to freezing onto warm 
grains processed through the J-shock front 
and traveling downstream along the mag- 
netic field lines. 

5. the asymmetry in optical and NIR proper- 
ties among the two lobes of the outflow is 
probably not caused by different pre-shock 
densities or shock velocities, supporting the 
hypothesis of higher extinction values to- 
ward the red lobe. 

6. the total J-shock cooling is compatible with 
the molecular outflow being accelerated by 
the fast neutral wind detected in Hi and CO. 

7. the whole flow area appears to be associ- 
ated with a faint PDR illuminated by BD 
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+30° 549, the source responsible for the il- 
lumination of the whole NGC 1333 nebula. 
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